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Two-jet event shape distributions, traditionally studied in the language of perturba- 
tive QCD, can be described naturally in soft-collinear effective theory. In this language, 
we demonstrate factorization of event shape distributions into perturbatively-calculablc 
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1. Introduction 

Hadronic jets produced in e + e~ annihilations provide a rich laboratory for probing 
both the perturbative and nonperturbative behavior of strong interactions described 
by Quantum Chromodynanhcs (QCD). Studying these leptonic collisions allows 
us to focus on strong effects only in the final state. Asymptotic freedom makes 
possible the reliable description in perturbation theory of processes at high center- 
of-mass energies Q, such as the prediction of dominance of two- and three-jet-like 
states in the hadronic final states in e + e~ annihilation and their respective cross- 
sectionsP The nonperturbative process of hadronization at the low energy scale 
Aqcd, however, introduces corrections which we do not yet know how to calculate. 
For some observables, such as the total hadronic cross-section, these corrections 
are small, of order (Aqcd/Q) 4 -^Ih less inclusive observables, these corrections are 
larger, on the one hand increasing uncertainty due to nonperturbative corrections, 
but, on the other, giving us an opportunity to measure their size and probe some 
of their properties even if we cannot calculate them exactly. 

One such class of observables is that of event shapes, numbers that depend on the 
distrbution of momenta of final state hadrons, probing in greater detail the structure 
of the jets than does the total hadronic or n-jet cross section. The study of e + e~ 
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event shapes has a long history in the framework of perturbative QCD (pQCD)P A 
two-jet event shape e is defined to range between and 1, with exactly back-to-back 
pencil-like jets corresponding to e = 0. Predicting the distribution of events near 
the region e = requires both perturbative resummation as well as an accounting 
of nonperturbative power corrections, which enter as powers of Aqcd / (eQ) due to 
nonpcrturbative soft or collinear radiation from jets. 

Many approaches have led to fruitful analysis of these power corrections, mostly 
based on the behavior of perturbation theory predictions for event shape distribu- 
tions. Renormalon ambiguities of order Aqcd/Q in the perturbative series for jet 
observable s led to the prediction of compensating Aqcd/Q nonperturbative power 
corrections! 4 l ^ | 7 | 8 | 9 | 10 | ll | 12 | 13 | 14 | These studies led to models of a universal ef- 
fective or dispersive infrared coupling replacing the strong couplin g at low scales 
implying universal power corrections to differe nt event s hapes! I land later to the 
model of dressed gluon exponentiation (DGE) P^^^ combining the resummation 
of perturbative logarithms and nonperturbative corrections from renormalons. Sec- 
tion 2 introduces the event shapes we choose to study and reviews the predictions 
of universality. 

Another powerful tool in the study of power corrections is f actorization, sepa- 
rating perturbative and nonperturbative effects in jet observables L12118119I20J j^vent 
shape distributions can be written as convolutions of perturbatively-calculable and 
nonperturbative functions, written schematically: 

= H ® J x ® J 2 ® S, (1) 

de 

where H describes the hard scattering, J\_2 describe the dynamics within the jets, 
and S describes soft radiation away from the jets. The soft function S can be viewed 
as a resummation of all powers of AQCD/(eQ) in the region e ~ Aqcd/Q, but we 
can also look over a larger range of Ae, and extract information about the low-order 
terms in AQCD/(eQ) from S. 

The separation of scales (hard, jet, and soft) that allows a factorization of the 
form of Eq. fl} suggests the formulation and use of an effective field theory descrip- 
tio n in c alculating the event shape distributions. The soft- collinear effective the- 
or jj2Tl22] (gCET) is the theory appropriate for this task. Originally applied mostly 
to problems in B physics, SCET also provides a natural description for the en- 
vironment of jet physics, an d has been used not only to analyze nonpe rturba tive 
power correction d^ 24 ^^ ! but also the improvement of event generator J22EH1 anc j 
the perturbative calculation of jet cross sections^! In this review we focus on the 
identification of the nonperturbative power corrections, and not so much on the 
perturbative issues of matching, running, and resummation. Section 3 provides a 
brief overview of aspects of the SCET formalism that we need in this discussion. 

In Sec. 4, we use SC ET to reproduce the factorization theorem £[]) for event 
shape distributions^^], already well-known from arguments 

in p Q Cr j l2 | 20 | 3 0| 

and 

identify more precisely what is meant by the schematic notation in (fl} . Wc will not 
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be able to provide a detailed description of the pQCD treatment of event shapes in 
this review, focusing instead on the SCET formulation. However, a comparison of the 
treatments in SCET and full QCD has not only elucidated tight connections between 
the two approaches, but has also inspired proofs of the universal relations among 
power corrections to different event s hapes which go beyond some of the assumptions 
implicit in the pQCD-based models! I I We review these developments in Sec. 5, 
and conclude in Sec. 6. 



2. Event Shapes 

2.1. Two- Jet Event Shapes 

The event shapes we shall consider are defined as functions of the momenta pi of 
the particles in the hadronic final state N in the collisions e + e _ — > N. They are 
"traditional" event shapes such as the thrust: 

T=-max^|t-p i |, (2) 



ieN 



the jet broadening: 



and the C parameter: 



B = ^J2\ ix p^ ( 3 ) 



ieN 



3 

a 

ieN 

and the more recently introduced angularities^^. 

T a = i^^sin Q 0i(l -cos 6 l ) 1 - a , (5) 
^ ieN 

where the angles 9i are measured with respect to the thrust axis. The angularity r a 
is an infrared-safe observable for a < 2, although our analysis of nonperturbative 
power corrections below will be valid only for a < 1. Between a = and a = 1, 
the angularities interpolate between the thrust (tq — 1 — T) and the jet broadening 
(ti = B). Adjusting the continuous parameter a to move through this family of 
similar event shapes will allow us to extract more information about nonperturbative 
physics from the power corrections than would be possible by studying individual 
event shapes separately. 

These event shapes can be rewritten in a common form, using the rapidities rji 
of final state particles, where rji = | In ^'^i ^ , z being along the thrust axis. For 
massless particles, rji = lncot(#i/2). The event shapes above take the general form: 

G = 7)£|P^|/e^)> ( 6 ) 

^ ieN 
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where the appropriate choice of the function f e yields the various event shapes. 
For the angularities, f Ta (T a ) = exp[— |?7,| (1 — a)], and for the C-parameter, fc = 
3/ cosh r^. The form ([6]) will be useful for studying the behavior of soft gluon radi- 
ation under Lorentz boosts along the thrust axis. 

2.2. Power Corrections from Perturbative QCD 

Analyses based on perturbative QCD predict a simple behavior for the leading 
power corrections to event shape distributions due to soft radiation. Many of these 
are based on finding the shift in the resummed perturbative event shape distribu- 
tions due to the (exponentiated) emission of a single gluon with small transverse 
momentum, k 1 - ~ Aqcd ■ The consistent prediction of such approaches is a universal 
power correction shift ing the perturbative predictions of the mean valued^ or even 

IQOI 

the entire distribution ^ of event shapes: 



where c e is an observable-dependent but exactly-calculable coefficient and A is an 
unknown but universal nonperturbative parameter. This universal structure persists 
up to 0{a 2 s ) in soft gluon emission, with the single-gluon shift, c e A/Q, just multi- 

IOOI 

plied by another universal number, the Milan factor A'ff^In approaches modeling 
the strong coupling with an effective infrared^ or dispersive coupling23 m the 
low-momentum regime, A depends on the IR coupling <5o(^/), where fii is the 
infrared-cutoff scale dividing perturbative and nonperturbative regimes, and the 
perturbative strong coupling a s (/x) at a higher scale. Fits to the shift in mean val- 
ues of various event shapes observed in data yield a value of a s (Mz) consistent with 
other extractions and a universal value of ao(fJ>i) at fXi = 2 GeV of about 0.5oP 
Good fits to the event shape distributions over the whole two- jet reg ion and a wide 
range of Q 2 have been obtained using the DGE cross-sections! I I 

Similar arguments applied to the angularity distributions predict a remarkably 
simple scaling behavi or for the entire soft function appearing in the factorized an- 
gularity distributions IMI 35 | 36 | distribution in r a takes the factorized form: 



where o~°j is the perturbatively-calculated r a distribution, and S a (£) is the soft 
function. Analyzing the contribution of low-momentum gluon radiation in the NLL 
resummed cross-section, the soft function is shown to exponentiate, 



transformed to moment space by S a {v) = J^ d^e v ^S a (^). According to this pre- 
diction, the coefficients A„ are independent of a, and the various angularity shape 






(8) 
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functions are related to the thrust shape function by: 

S a (v) = S (v)^ . (10) 

This relation has been tested using simulated data from PYTHLAp^, and awaits 
tests against data from LEP. 

The predictions of these perturbative QCD-based models for the soft power 
corrections are quite successful, but as they rely on the behavior of perturbation 
theory for low-order gluon emission, seem not to capture fully the nonperturbative 
effects to all orders of soft gluon emission. To move beyond this limitation, we will 
find useful an approach that clearly separates the perturbative and nonperturbative 
physics, provided for us by the tools of effective field theory. 

3. Soft-Collinear Effective Theory 

In this section we provide a brief but mostly self-contained overview of the formalism 
of soft-collincar effective theory. 

3.1. Soft and Collinear Modes 

In a two-jet event, there are several distinct momentum scales: the total center-of- 
mass energy Q, the scale of nonperturbative radiation Aqcd, and the intermediate 
scale ^QAqcd describing the typical transverse momentum of particles in the jets. 
The existence of this hierarchy of scales leads to the construction of an effective 
field theory, th e soft - collinear effective theory as an expansion in the parameter 
A = \/Aqcd/Q !^ The modes in this effective theory characterized by the scal- 
ing of the light-cone components of their momenta, p = (n-p, n-p^p^), with respect 
to back-to-back lightlike vectors n, n. Collinear momenta scale as p c ~ Q(A 2 , 1, A), 
and ultrasoft (usoft) momenta scale as p us ~ Q(A 2 , A 2 , A 2 ). The theory with this 
collinear scaling is called SCETi, and is appropriate to describe jets whose con- 
stituents have typical transverse momenta of the order p-^ ~ ^/QAqcd- Narrower 
jets may be des crib ed in a theory SCETn in which collinear momenta scale as 
Pc ~ <3(A 4 , 1, A 2 )P^ We will study event shapes dominated by jets wide enough to 
be treated in SCETi, and the shortened abbreviation "SCET" will always refer to 
SCETi in the remainder of this review. 

The fields in the SCET Lagrangian are collinear quarks and gluons, £ n and A n , 
and ultrasoft quarks and gluons, q us and A us . For the collinear modes, we split their 
momenta into two parts, a large label piece and a residual piece: p n — p n + k, where 

P%=n-PnY+Pn (11) 

contains only the 0(Q) piece of n ■ p n and the 0(QX) piece of leaving the 
residual momentum k of order QA 2 in all components. The collinear and soft fields 
are defined by factoring out a phase removing the dependence on the label momenta 
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from the full theory fields q and A: 

q(x)=^2e- i P- x q n>p (x) + q as (x), A(x) = ^ e'^* A n>p (x) + A us (x). (12) 

The fields q n ,p(x) and A n ^ p (x) now have fluctuations in spacetime described by 
momenta of order QX 2 = Aqcd about the corresponding label momentum p. The 
restriction of the sums to nonzero labels is crucial, to avoi d d ouble-counting the 
momenta already covered by the usoft modes q us and A„ S P3 It is useful to de- 
fine operators^ which pick out only the label components of the momentum of a 
collinear field </>„ )P : 

Vp,4> n ,p = P l i4>n,p, 'Pfi^p = ~P^n,p- ( 13 ) 

We will also use the shorthand P = n ■ V, and V\. 



3.2. The SCET Lagrangian 

To complete the construction of the SCET Lagrangian, the components 

t _ # „ _$4_ nA s 

are projected out of the field q n ,p- Writing the QCD Lagrangian in terms of these 
fields, it can be shown that £„ jP is massless while S„ p has an effective mass n ■ p ~ 
0(Q), which is thus integrated out of the theory. Keeping only the terms of leading 
power in A, we obtain the leading-order Lagrangian of SCET: 

^SCET = £cq + £cg + £-us , (15) 

where the collinear quark Lagrangian L cq ia^l 



in ■ D us + gn ■ A n , q + (fj_ + g4n,a)W n ^W£(T ± + g4i Q 



V 



(16) 

where W n (x) = Pcxp[ig f_ dsn ■ A n (ns + x)] is a Wilson line of collinear gluons; 



the collinear gluon Lagrangian ig®2 

Ccg = 2^ Tl ' ^ + gK ^ ^ + + ^ ( 17 ) 

where iV 1 = + 7 -^-in- D us and C 9 C ^ contains the collinear ghost and gauge-fixing 
terms; and the ultrasoft Lagrangian, 

C us = q us ip us q us - \ Tr G^G^ + C 9 J- , (18) 

is simply that of QCD. The ultrasoft covariant derivative iD£ s = id^ + gA£ s con- 
tains only the ultrasoft gauge field, the ultrasoft gauge field strength is G^ v — 
(i/g)[D£ s ,D" ls ], and contains the ultrasoft gauge fixing and ghost terms. In 

the collinear Lagrangians, the label momentum phase factors and sums over labels 
are implicit. The se simply have the effect of enforcing label momentum conservation 
in each termPSl 



I 



I 
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3.3. Decoupling Soft and Collinear Modes 

In the leading-order SCET Lagrangian, usoft quarks do not interact with collinear 
modes, and the coupling of usoft gluons to collinear modes is particularly simple, 
with only the n ■ A us component appearing. These interactions can be ma de t o 
"disappear" from the collinear Lagrangians by the following field redefinitions!^ 

^=4C- A n =Y}A' n Y n , (19) 

using the outgoing Wilson line of ultrasoft gluons, 



Y n (z) = Pexp 



ig / dsn ■ A us (ns + z) 



(20) 



The redefinition |T9]) also implies W n = Y\W n Y n . Since in ■ D US Y^ = 0, making 
these substitutions for the collinear fields in C cq and C cg removes their couplings 
with usoft gluons: 

£ n ,p>in ■ D us Cn^p = £ n yin ■ d£' np) (21) 
Tr { \iV» + gA* t9 , iD v + gA» tq ,} 2 } = Tr { + gA% q , iVfa + 5 < g ,] 2 } , (22) 

where i'Dyy* — 7-*^ + if*™ ' This absence of interactions between the collinear 
and ultrasoft sectors simplifies proofs of factorization at leading order in the SCET 
expansion. Effective theory operators built out of collinear fields will always have 
ultrasoft Wilson lines accompanying them. 

4. Factorization of Two-Jet Event Shape Distributions 
4.1. Matching the Two- Jet Current 

The two-jet event shape distributions in e + e~ — ► 7* — > N are calculated from 
matrix elements of the current = qj^q in QCD: 

Te = W ^ m jM(0) |0) L ^ (2 ^ )4<54(0 ~ PN)6( - e ~ e(iV)) ' (23) 

where L M is the leptonic part of the amplitude. We want this current to create to 
create states with two hadronic jets back-to-back in the center-of-mass frame. So, 
we will match this current in QCD onto operators in an SCET for collinear fields 
in two separate light-like directions^l, m and ri2- 

Jni,n 2 ~ Xni,p 1 J±Xn 2 ,P2 1 (24) 

where we have introduced the jet field x-n = W^ n . The presence of the collinear 
Wilson line is required by gauge invariancePSl The operators satisfy the matching 
condition: 

(i^)QCD(M) = C ni „ 2 (pi ■p2;M)(Jn 1 ,n 2 )sCET(M), (25) 

a Thc field redefinition may be performed with incoming Wilson lines, with s = —00 to 0, depending 
on whet her the collinear particles to be described are incoming or outgoing in the hard scattering 
or decay.HH 
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where C is the matching coeffcient, dependent on the matching scale /i. A sum over 
ni,2 is implicit. The matrix elements on both sides are taken between the same 
states, calculated in QCD and SCET, respectively. Here, our currents must produce 
a back-to-back qq pair in the final st ate, which requires ri\ = fi2- To 0{a s ), the 
necessary matching coefficient iJSHEE] 



47T 



n i 1 2 

Y + ln 



3 In 



(26) 



P\-PiJ \Pi-P2, 
Qn/2 andp2 = —Qn/2 (t; n ,p 2 creates an antiquark 



so pi ■ p2 = —Q 2 - (For the current in deep inel astic 

Q 2 .) 4 3 The 



The collinear labels must be 
with label momentum 

scattering in the Breit frame, both labels are positive, and pi ■ p2 
matching coefficient simplifies at the scale u = Q: 

a s (Q)C F f q 7ir 2 

The matching coefficient at another scale fj, can be obtained via r enorm alization 
group evolution. The anomalous dimension of the two-jet operator iJSZEH] 



C nfi (-Q 2 ;Q) = 1 



3l7T 



(27) 



72 (m) 



log 



(28) 



The accuracy of the effective theory prediction for the event sha pe distrib ution can 
be improved by including operators with higher numbers of jet cP^^^ (collinear 
fields in more directions). 



4.2. Factorizing the Distribution 

With this matching, the SCET expression for the event shape distribution becomes 



da 
de 



KM 

2Q 2 



J2 I (N\ TXn, Q Yn7lY£xk,-Q (0) |0) L M (2^) 4 <5 4 (Q -p N )6(e ~ e(N)) 



N 



(29) 

where we have made use of the collinear fields x'n = ^n^n redefined according to 
Eq. (fi"9|) . so that the soft gluons appear in Wilson lines in the operator. The time- 
ordering which was implicit in the matrix element in Eq. (|23|) is written explicitly 
in Eq. ([29| since the Wilson lines contain fields at many different times. 

To factorize the matrix element in (|29p into purely collinear and soft parts, we 
can split t he final state \N) into collinear and soft sectors, \N) = \J n ) ® \Jn) ® 
\X U ) 12413H Particles in the final state are assigned to these sectors according to 
their label momenta — particles with zero label momentum belong to the sector 
.X^PHl Since there are no interactions in the Lagrangian between the collinear fields 



,, fl , A' n fi and usoft gluons A US7 the matrix element in (|29|) factorizes to give: 



^ = §2 E K^^JTx'^x^-qWIO)^! 2 ^^^^^^^^)]^} 

x (2n) 4 d 4 (Q - PJn -p Jn - Px J6(e - e(J„) - e(J n ) - e{X u )), 

(30) 



I 



l 
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where the trace is over colors. We can write this factorized cross-section as a con- 
volution of purely collinear and soft functions: 



— = dejaj(e.j)S e (e- ej), (31) 



where aj is the collinear cross-section 

CI 2 2 
CTj ( ej ) = 2Q2 ^ \( J n J n\Xn,Ql±X'n^Q(°) \°) L p.\ 



x (2tt) 4 5 4 (Q -p Jn -pjJSiej - e(J n ) - e(J n )), (32) 



which is equivalent to the perturbatively-calculated distribution do~/de\pT, and S e 
is the soft function 

S e (e) = -Lti-]T \(X u \T[Y n Y£] \0)\ 2 6(e - e(X u )). (33) 

To separate the two functions we neglected t he c ontribution of the usoft momentum 
in the momentum-conserving delta functiorl^U in Eq. (I32p . which is a subleading 
contribution in Aqcd/Q, and took the usoft momentum in the soft function (|33() to 
be unrestricted. This is correct because the high momentum cutoff of the usoft in- 
tegrals can be viewed as parametrically larger (in 1/A ) than typical usoft momenta, 

1041 

so the upper limits are infinite in the limit A — > 0!^ 



4.3. Eliminating Double- Counting 

This procedure raises the question of double-counting. The soft function contains 
radiation emitted along the jets, overlapping with collinear radiation whose en- 
ergy is small. This double-counting is avoided in SCET by always implementing 
the restriction to nonzero label momenta in the sums in (|12p for collinear modes 
when calculating collinear phase space or loop integrals. This procedure is kno wn as 
zero-bin subtraction and has important consequences in calculations in SCETpSl A 
related procedure in factorization theorems of perturbative QCD is the introduction 

of eikonal jet functions and eikonal cross- sections which are subtracted out of the 

101 1 

naively-defined jet and soft functions to eliminate double-counting. Interpreting 
these functions in the language of SCET leads to a simple demonstration of the 
equivalence of th e tw o su btra ction methods, shown to leading order in A and all or- 
ders in a s in 

Ref M Re ,Hl also illustrated this equivalence by explicit calculation 
at fixed orders in a s . 



4.4. Running the Jet and Soft Functions 

The collinear and jet functions in Eq. (|31[) depend on scale. As we saw in Eq. (|26p . 
the two-jet matching coefficient C„ fi (/i) contains logarithms of /i/Q, and so is nat- 
urally evaluated at the scale a = Q. The collinear jet function and soft functions 
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are naturally evaluated at scales /i c ~ QX and [i s ~ QX 2 . The running of the func- 
tions among these different scales requires knowing the anomalous dimensions of 
the collinear and soft matrix elements appearing in the jet and soft functions!^ 
For some event shapes, nonperturbative contributions to the collinear jet function 
at the scale fi c ~ QX 2 become important. In this case, we need to know how to run 
the collinear matrix element down to this scale. Between the scales QX and QX 2 , 
the scaling of the collinear momenta changes from that of SCETi to that of SCETn, 
and we must match the first theory onto the second. The fields and Lagrangian of 
SCETn remain the same, only the scaling of the collinear momenta changes, and 
collinear matrix elements in SCETn are nonperturbative. 

We can estimate the size of nonperturbative power corrections in the collinear 
and soft functions for a particular event shape by considering the contribution of par- 
ticles with transverse momenta of order p- 1 ~ Aqcd to these functions in SCETn- 
Rapidities of collinear particles in SCETn are of order \rj\ ~ i |ln(n -p/h ■ p)\ ~ 
ln(l/A 2 ), and those of soft particles are of order r\ ~ 0. In the case of angu- 
larities, the factor / Ta (r a ) = e I^K 1 **) for collinear particles is then of order 
^2(1 -o) _ (Aqcd/Q) 1- ", and for soft particles, order 1. Then the nonperturbative 
collinear radiation contributes an amount of order |pj_| e - ' 1 ^ 1-0 ) ~ (Aqcd/Q) 2- ", 
and the soft radiation an amount of order Aqcd/Q, to the event shape T a . Counting 
more carefully, we find the collinear and soft power corrections to the r a distribu- 
tions to be in powers of (AQCD/(iaQ)) 2 ~ a and AqcdAtciQ), respectively. Then for 
a < 1, the soft power corrections are dominant, and we can ignore the collinear 
power corrections to first approximation. For a > 1, the collinear power corrections 
dominate. For the broadening B they are of the same order. Furthermore, in this 
region, there are also other sour ces of po wer corrections we have chosen to ignore, 
such as recoil of the thrust axisP^ESESI p or these reasons we limit our treatment 
to r a with a < 1 as well as other event shapes such as the C-parameter for which 
the soft power corrections dominate. 

5. Soft Power Corrections 

We now turn our attention solely to the power corrections coming from the soft 
function in the factorized event shape distribution (|3ip . By studying the energy flow 
and boost properties of the soft radiation, we will be able to show the universality 
of the power corrections to the mean values of event shapes. This universality will 
be violated in shifts of higher moments, but for the angularity distributions, we will 
be able to extract useful information about the nature of the soft radiation from 
the details of these violations. 

5.1. Energy Flow 

The delta functions of the event shape variables in the jet and soft functions (f3"2"|) 
and (|33p can be expressed in terms of the transverse momentum flow operator, 
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£t(v)i defined by its action on a state \N): 

£ T (r 1 )\N) = Y / \pt\S(v-V l )\N), (34) 

iejv 

adding up the magnitudes of transverse momenta of particles with rapidity rj. These 
operators are simply related to the energy flow operators^ £{h) 1 which act on states 
according to: 

£(n)\N) = Y,E t 6 2 (h~n l )\N), (35) 

measuring energy flowing in the direction of the unit vector h. With S 2 (h — hi) = 
5(cos6 — cos8i)S((f) — 4>i), and rj = In cot (0/2), the two operators are related by: 

£t(v) = —^r- I d q b£{h). (36) 
cosh rj Jo 

The energy flow operator c an be repre sented in terms of the energy-momentum 
tensor on a spher e at infinity-^^^IS] anc j j s closely related to the energy-energy 
correlation J 52 ! 5 ^ !, 

Using the operator &r(r]), the soft function can be rewritten: 



1 - / 1 f°° 

S e (e) = — TrJ2(0\ f{YnY^}d [e-- J_Jr,£j 



( V )f e ( V ) )\X u )(X u \T[Y n Y 7 t}\0) 



(37) 

Since there is no restriction on the usoft states in the soft function S e (e) in Eq. I|37[). 
we can perform the sum over states to obtain 

S e (e) = ^- Tr (0| Y^S d V £ T (v)fM) Y n Y fl |0) . (38) 

The time- and anti-time-ordering operators have been removed by using Wilson 
lines Y nyn in the anti- fundamental representation of SU(N), and the space- like 
separation (ns — hs') 2 — —Ass' < of the n and h fields!^ 



5.2. Boost Invariance and Universality 

Inserting the soft function ([38| back into the convolution ([3l"j) , we can take moments 
of the full distribution. The first moment (e)px of the perturbative distribution 
cj(ej) is shifted by the first moment of the soft function: 

(e) = (e) PT + £(e) N p, (39) 

where 

£(e) N p = q J dr, U(ri)— Tr (0| yIy^t^Y n |0) . (40) 

We can make use of boost properties of the operators in the soft matrix element 
to simplify this quantity. Insert factors of 1 = U^U^, where is the operator 
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implementing a Lorentz boost by a rapidity rf along the thrust axis, around each 
operator in the matrix element: 

(0| U/U n Ylu-, l U^U-^S T {r,)U^U^ |0) . (41) 

Now we use the transformation properties of each object in this matrix element: 

Efy|0) = |0), U n ,Y n Uj 1 = Y n , U v ,£ T {'n)U-, 1 = £T{ri + i 1 l ), (42) 

where the boost invariance of Y n can be seen from: 



Y n = P exp 



ds n ■ A us (ns) 



— » Fexp 


ig \ 




. Jo 



dse v n ■ A us (e n ns) 



= Y„ 



(43) 

by a simple change of variables s' — e~ v s to give the final equality. This is a mani- 
festation of the reparameterization invariance of SCETJ^The boost transformation 
of £t(v) is seen easily from its definition (f3~4"|) : 

U^&r^U-, 1 \N({ Vl })) = U^&riri) \N({ Vt - ?/})} 

= J2 u, r \pi\S( V - m + rf) \N({ m - V})} 

i£N 

= £ T (v + v')\N({Vi}))- (44) 

This means that the soft matrix element in Eq. (|40j) is independent of the rapidity 
chosen as the argument of £t(v) : 

(0| YlY^£ T { v )Y n Yn |0) = (0| Y j n Yj£ T (r, + V ')Y n Y n |0) = A (45) 

This matrix element does not depend on 77 at all, and so is just a constant A. We 
thus take it out of the integral in Eq. (|40|) . leaving 



where 



A 

(e) = (e)p T + c e — , 



drif e (r,), A={£ T (0)), 



(46) 



(47) 



the Wilson lines and color trace being understood in the expectation value of £>r(0). 
We have chosen the argument of £t(ji) to be 0, but could have chosen any rapidity. 
For the angularities, c Ta — 2/(1 — a), and for the C-parameter, cc = 37r. 

These observations reproduce the prediction ((7|) for the universal shift in the 
mean value of event shapes. We have not made any assumptions about single or 
low-order gluon emission in the final state. This result relied rather on the feature 
of the SCET Lagrangian at leading-order in A that only the n ■ A us component 
of usoft gluons couples to ri-collinear particles (n ■ A us to n-collinear particles), so 
that these couplings can be removed from the Lagrangian by the redefinition (|19p 
of the collinear fields, and on the boost invariance of the usoft Wilson lines. The 
observation that the e-dependent effects of soft radiation enter only through the 
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integral in (|47|) for c e was already noted in the derivation of the Milan factorPSl 
Here, we have obtained in addition a field-theoretic interpretation for the quantity 
A, and, at leading order in the SCET expansion in A, demonstrated its universality 
exactly to all orders in soft gluon emission. 

Through same line of reasoning we discover that the prediction of a universal 
shift in the full distribution is spoiled by nontrivial correlations in the energy flow 
of soft radiation. For example, the shift in the second moment of the perturbative 
distribution due to soft radiation is given by: 

£(e 2 >NP = / dr H / d m feMfeim) (Q| YlY^£ T ( m )£ T ( V2 )Y n Y n |0) . (48) 

V J — oo J — oo 

Boost properties of the soft operators can be used to eliminate dependence on one of 
the rapidity variables but not both. Thus this shift does not simplify to (c e A/Q) 2 , 
as a universal shift ([7]) of the full distribution would require. 

5.3. Power Corrections to Angularities 

More information can be gleaned from the behavior of power corrections to the 
angularity distributions. First, we transform the soft function S e (e) to moment 
space via the Laplace transform, 



SM=J dee- ve S e (e) = ^exp 
which exponentiates in terms of cumulant moments 



v 

QJ-c 



dr/£T(v)fe(v) 



S e (^) = S e (0)exp 



|M-i)"((/>^- w 



(49) 



(50) 



where the cumulants are defined according to ((X)} = (X), {{X 2 }} = (X 2 ) — (X) 2 , 
and so on. 

For the angularity shape function, we may compare this form to the prediction 
(|9j) based on the behavior of resummed perturbation theory. The parameters A„ in 
Eq. ([9]) correspond here to the cumulants of energy flow operators: 

j^Ua) = {(f"j*frto>~ Ml ~ a) ) n }> ( 51 ) 

where we now allow a possible a dependence in the coefficients A„ . For n = 1, we 
find simply 

1 2 

■ -Ai = - A, (52) 



1 — a 1 — a 

and so Ai is truly a- independent. However, beginning with n = 2, the coefficients 
A n depend, in general, on a. For example, 

-A 2 (a) = -/ dm dme-^- a ^ + ^{{&r(m)£r{m)))- (53) 



1 - a v ' 2 



X 
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We again insert Lorentz boosts into the matrix element, and this time are able to 
eliminate dependence on one but not both rapidities, and perform one integration: 
1 1 r°° 

" -^(a) = ^ \ / *?[! + (!- <0 \v\}e- M(1 - a) ((£T(0)£T(v))), (54) 



I- a w 2(1 - a) 

so correlator still depends on the rapidity gap between the two insertions of £t- 

If the two-point correlations of energy flow ((£t(0)£t(v))) are negligible for ra- 
pidities 1 77 1 3> jzr, then the coefficient of the matrix element in (|5"4"|) can be expanded 
about r\ = 0, 

A 2 (a) - i |~ drj |l - - a) V } 2 + ±[(1 - a) \ V \f + ■ ■ ■ } <(£ T (0)£ T (r?)}}. (55) 

Thus the a- independence of A2 (a) is tied to the assumption of negligible correlations 

1041 

in energy flow of soft radiation!^ 1 

The violation of the scaling rule (flC)]) can itself be used to deduce information 
about the correlations in the final state soft radiation. The quantity, 

1 r°° 

C 2 (a) = - d* ? e-l"l( 1 - B )((fir(0)fir(» 7 ))), (56) 



is essentially a Laplace transform of the two-point correlator to the variable a. 
Equation (|54|) relates 02(a) to the quantity \2(a)' 

C 2 (a) - dln{ i_ a f 2(a) - A a (o), (57) 

whose solution is 

C 2 (a) = (1 - a) ^-^X 2 (a% (58) 

assuming the correlator vanishes for a — > -co. Thus, the observation of a- 
dependence in the coefficient A2 (a) in the expansion of the angularity shape function 
([9]) yields direct information on correlations of the energy flow of soft radiation at 
separated rapidities. This information can, in principle, be extracted from data in 
e + e~ collisions at LEP. 



6. Conclusions 

Event shapes provide a useful laboratory for the study of perturbative and non- 
perturbative properties of hadronic jets. The continuous adjustable parameter a 
in the angularities provides an extra source of information about the nature of 
the hadronization process in jet formation. Soft-collinear effective theory provides 
an organized framework in which to understand factorization, running, and power 
corrections that have traditionally been described in the language of full QCD. In 
particular, the decoupling of collinear and usoft modes in the leading-order La- 
grangian facilitates factorization, and the form of the soft function and its boost 
properties lead to a proof of the universality of power corrections to the mean values 
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of event shapes, as well as the extraction of useful information from the violation of 
universality or scaling of power corrections to higher moments. SCET systematizes 
the approximations that lead to useful predictions, and sets up the framework that 
would allow one to calculate corrections, such as those suppressed by A or that 
come from including operators with more jets in the matching from QCD. A fuller 
understanding of power corrections to event shapes in e + e _ collisions helps point 
the way to the eventual understanding of jets produced in hadron collisions and the 
nonperturbative dynamics of QCD. 
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